Creating Steady Atom-Molecule Entanglement and Coherent Molecular output by 
Demkov-Kunike Type Non-adiabatic Transition 



P. Zhang 1 - 2 , H. T. Quan 1 and C. P. Sun 1 - 3 ' a ' b 
1 Institute of Theoretical Physics, the Chinese Academy of Sciences, Beijing, 100080, China 
1 Department of Physics, the University of Hong Kong, Hong Kong, China 
3 Department of Physics, Nankai University, Tianjin, 300071, China 

We, based on the photoassociation of fermion atoms into bosonic molecules, propose a scheme 
to create the steady entanglement between the atom state and the molecule state inside an optical 
lattice. The stability of entanglement state is guaranteed by sweeping the frequency of Ramman laser 
beam through resonance according to the second Demkov-Kunike (DK2) nonadiabatic transition 
model [1,2]. The probability amplitude of each components can be precisely controlled by adjusting 
the sweeping parameter. Considering the loss of molecule, the steady coherent molecular output 
can also be obtained in each lattice site. 
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In recent years, the photoassociation of fermionic 
atoms forming bosonic molecules [3-6] has aroused great 
interest. In this process, two atoms in condensate are 
illuminated by two Raman laser beams. This atom pair 
absorbs an photon from one laser beam and emits a pho- 
ton to another beam. As a result of stimulated Raman 
effect, there is a transition from the two atom state to a 
single molecule state. Most recently, there are many ex- 
periments to demonstrate this photoassociation process 
and the Bose-Einstein Condenstation (BEC) of molecular 
has been observed [5] . There also exists another physical 
mechanism, Feshbach Resonance to create the molecu- 
lar from the very cold atoms even in situation of BEC. 
When the external magnetic held of an atomic system 
is changed, an entrance interaction channel and a closed 
channel can be resonantly coupled and the population 
of a two atom bound state, i.e. the molecular sate, can 
be resonantly enhanced. In recent experiments, the dis- 
appearance of atoms and the creation of molecules from 
BEC have been observed [7-9]. 

In both the two cases, an interesting question is 
whether one can controll exactly the rate of producing 
the composite particle-molecule and even obtain the co- 
herent output of molecular [10]. If the dissipation mech- 
anism is not considered in this photassociation process, 
there must exists a Rabi oscillation between atomic state 
and molecular state and then there is only a "pulse pro- 
file" of molecular productions. In this case, we can not 
obtain a steady population of atom and molecular state. 
This phenomenon is very similar to " pulse atomic laser" . 
In the experiment of atomic laser, an crucial technology 
to overcome this problem is to sweep the frequency of 
radiation field (r.f) [11] that couples the trapped BEC 
state to the un-trapped state. This physical mechanism 
avoiding the pure pulse profile due to Rabi oscillation 
can be understood in the rotation picture. Within this 
representation, the sweeping frequency will shift the orig- 
inal energy level of atom, Then the sweeping derives the 
system to cross the degenerate point of energy level in a 



non-adiabatic way. Then, the Landau-Zener approxima- 
tion [12] can be utilized to obtain the steady populations 
of BEC atoms in the different energy levels. In this arti- 
cle, we will invoke this idea of sweeping frequency to ob- 
tain a steady atom-molecule entanglement state, i.e. the 
coherence superposition of atomic and molecular num- 
ber states with steady population of molecules produced 
from the atoms in photoassociation. 

We notice that the Landau-Zener model for sweeping 
the frequency is precisely correct, only when the detun- 
ning in the effective two level Rabi model linearly varies 
with time t from — oo to +oo. Actually, the linear func- 
tion of time is only an approximate description of the 
variation of the realistic detunning, because in practical 
cases the experimental parameters are always adjusted 
as smooth functions of time which have finite values in 
the limit t — > ±oo. For instance, in ref. [13], the hyper- 
bolic tangent function is considered as a more practical 
description of the experimental behavior of the atomic 
Stark energy difference controlled by laser beams. The 
validity of the Landau Zener model is also analyzed nu- 
merically in the same reference. 

In this article, we adopt a feasible sweeping way de- 
picted by the second Demkov-Kunike (DK2) [1,2] model. 
In this model, the detuning has finite asymptotic values 
in the limit t — > ±oo. When we sweep the detuning in 
the DK2 way, the finite asymptotic values of detuning 
are experimentally implementable. 

In our study, the model Hamiltonian for the atom- 
molecule photoassociation inside a single site in a optical 
lattice can be written as [14]: 

H = n b tfb + -U h b% [6 f 6- l] (1) 
+Q,f<j z + XJ x b^bo z + xb cr_ + x&cr + . 

Here, b is the annihilation operator of molecule and x 
is proportional to the two-photon Rabi frequency asso- 
ciated with the two laser beams [14]. Without loss of 
generality, in the following discussions, we assume x to 
be real, fib and flf are defined as 
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n b = u b + a + u x , 

where Wf, is the atomic energy, u)f is the molecular energy 
and A the detunning of the two laser beams. The terms 
contains Uf, U x and U b denote the energies of two body 
interactions between atoms, atoms and molecules, and 
molecules. In the Hamiltonian (1), the Anderson's qusi- 
spin maps [15] defined by 



C\C 2 , + 



c\ci 



c\c 2 
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have been used. Here, Cj (i = 1, 2) is the annihilation op- 
erator of the atomic state It is easy to see that the op- 
erators <r± are defined as the qusi-spin flips <j + = l e ) (ffl; 
c- = Iff) (e| and <r z = \e) (e| - |g) (,g| where |g) = |0) a 
is the atom vacuum state and |e) = c+cj |0) o is the two 
atom state. In the mathematical formulism, the Hamil- 
tonian (1) is just a Jaynes-Cummings kind Hamiltonian 
with a nonlinear detuning [14]. The laser enhanced com- 
bination process depicted by the interaction Hamiltonian 
can be illustrated in Fig 1. 

FIG. 1. The photoassociation process of the atoms in opti- 
cal lattice. In every lattice, the molecule state \M) (the two 
atom bound state) is coupled with the atom state \F) (the 
two atom free state) via the two Raman beams u\ and LJ2- 

In order to use the DK2 model, we assume that the 
detunning A of the two Raman laser beams varies with 
time as a hyperbolic tangent function: 



A (t) = 2uj f - oj b + U f - 2fctanh 



(3) 



It is obviously that the asymptotic values of A (t) in the 
limit t — ±oo are determined by the parameter fc and the 
speed of the variation of A (t) depends on T. Therefore, 
in the limit t — ► ±oo, the A has finite asympotic values 
2uf -uj b + Uf ± 2k. Since |tanh(±2)| = 0.96 w 1, the 
practical operational time is about AT. 

Obviously, the sub- 

space expanded by {|e, n b ) , \g, n b + 1)} is an invariance 
subspace V n of H in Eq. (1). In other words, the Hamil- 
tonian can be diagonalized as the block-diagonal matrix 
and an arbitrary superposition of \e,n b ) and \g,n b + 1) 

in V n can only evolve into V n itself. Here, we define 
— - + 

\e,n b ) — |e) <g> (rife!) 2 6'" b |0) b as the number state cor- 
responding to two atoms and n b molecules. Here, |0) b 
is the molecule vacuum state. \g,n b + 1) has the similar 
definition. Then a block of H can be expressed as a 2 x 2 
matrix: 



H{n b ,t) 



a + k tanh 



where we have used Eq. (3). Here, a and c are func- 
tions of n b and are defined as a (n b ) — n b (U x — U b /2) , 
c(n b ) = xV^bTT. 

Since \ is proportional to the Rabi frequency which is 
slowly varying with the frequency of the Raman beams, it 
is reasonable to assum9. 2847701 x 10~ 24 JT^'e c{n b ) = 
xV n b + 1 to be a time- independent constant. The in- 
stantaneous eigenvectors of the matrix H (n b ,t) can be 
noted as \$+(n b ,t)) and |$_ (n b , t)). It is apparently 
that in the limit t — ► ±oo, the matrix H (n b ,t) has 
asymptotic expressions 

H (n b , ±oo) = [a (n b ) ± k] |e, n b ) (e, n b \ + c (n b ) |e, n b ) (g, n b + 1| 
- [a (n b ) ± fc] \g, n b + 1) (g, n b + 1| + c (n b ) \g, n b -I 

Therefore, the eigenvectors |$± (n b ,t)) also have asymp- 
totic expressions when t —> ±oo. We note them as 
|$± (n b , ±oo)). Generally speaking, \&±(n b ,±oo)) are 
superposition states of \e,n b ) and \g,n b + 1). It can be 
seen in the following discussion that we are only inter- 
ested in the special cases that |<1>± (n b ,±oo)) can be ap- 
proximated as \e,n b ) or \g, n b + l). Then the general 
expression of |$± {n b ,t)) is not needed here. 

The instantaneous state |\& (n b , t)) of the quantum sys- 
tem is the solution of the Schorcdingcr equation 
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and it can 
|*± (n b ,t)) 



H(n b ,t)\V{n b ,t)) 
be expanded as the superposition of 



|* (n 6 ,t)) = Di(n b ,t)\^i(n b ,t)). 
»=+,- 

Obviously, \D± (n b , t)\ 2 are the populations in the states 
|$± {n b ,t)) at any instantaneousness. In practical cases, 
the asymptotic value \D± (n b , +oo)| 2 in the limit t — > 
±oo is usually of interest. If the quantum system is ini- 
tially prepared in the state 

|* (n b , -co)) = |$+ (n 6 , -co)) , 

the values of \D± (+oo)| 2 can be obtained with direct 
calculations [1,2]: 



\D- (+oo)| 



2 _ sinh [irTE+] sinh [7rTS_] 
sinh [irTEa] sinh [ttT^] ' 



(5) 



i^+(+^)i 2 = 1-1^- (+<x>)r 



where the complicated notations i? a , i? e and E± are de- 
fined as functions of the physical parameters a, b and 
k : 



|e, n b ) (e, n b | + c |e, n 6 ) (g, n b + 1| 



£ a (a, fc, c) 



a, fc, c) 



(a - fc) 2 + c 2 



(6) 



(a + fc) 2 + c 2 



a + fc tanh [ — 



Iff, ^fc + 1) (ff, m + 1| + c n 6 + 1) (eE^{ a , fc, c) = fc ± 



E P — En 
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In the above calculations, we have refered to the main 
result in the DK2 model. 

To create the atom-molecule entanglement state from 
the direct product of the atomic and molecular number 
sate \e,n b ), we assume the amplitude k is large enough 
so that the conditions k ± a {rib) » c (rib) are satisfied. 
Since a is proportional to rib and c is proportional to 
.y/nh, this condition implies that rib is small enough. In 
this case, when t = ±00, since k ± a (rib) >> c(rib), the 
diagonal terms of H (rib, ±00) is much larger than the off 
diagonal terms. The eigenstates of H (rib, ±00) can be 
approximated as 

|$+ (n b , -00)) ~ |$_ (n b , +00)) ~ \e, n b ) 
|$_ (rib, -00)) ~ |$+ (n b , +00)) ~ \g, n b + 1) 

This is an analogy of large detunning case when a two 
level atom interact with a classical radiation field. There- 
fore, in the limit t — > ±00, there is almost no atom- 
molecule transition and the mode of the probability am- 
plitudes of the states |e, n b ) and \g, n b + 1) do not change 
with time. This phenomenon has been shown in the ex- 
periment of photoassociation [5] . The mixing of the two 
states mainly happens in the nearby of the "resonance 
point" where \a + fctanh (t/T)\ < c. 

In this case, if the system is initially prepared in the 
state \e,nb) when t = —00 and A is changed with time 
according to condition (3), the system will evolve into 
the superposition state 



Di (n b ) \e, n b ) + D 2 (n b ) \g, n b + 1) 



(7) 



when t = +00. Apparently, this quantum state can be 
considered as a atom-molecule entanglement state. It 
follows from Eq. (5) that 



\Di (n b )| 2 



sinh [ttTE + (a, k, c)] sinh [irTE- (a, k, c)] 
sinh [irTE a (a, k, c)] sinh \ixTE e (a, k, c)] 



(8) 



where the functions E a (a,k,c), E e (a,k,c) and 
E± (a,k,c) arc defined in Eq. (6). The probability of 
creating one molecule from two atoms can be expressed 
as P = 1 — |Z?i| 2 . Especially when nb = 0, P describes 
the molecular vacuum production. 

It is easy to prove that when k is large enough so that 
e Tk » 1, l-Di (n b )\ 2 in Eq. (8) can be approximated as 



|L>i(n 6 )| «exp 



exp 



Tc 2 



(9) 



-2tt 



![2fctanh£] l*= 



to 



Here, t which satisfies a + k tanh (to/T) = is the cross 
point of the diagonal elements of the Hamiltonian in Eq. 



(4). Therefore, the approximated expression of \D\ (nb)\ 
in Eq. (9) is just the result given by Landau Zener for- 
mula. 

Now we are in a position to analyse the influence of 
the parameters T, nb, k, U x and Ub on the molecular 
probability P. It is convenient to express the probability 



_ sinh[7rTg_ 



P as P = l — fg where / = 
Using the simple relationships 



iYiWTE a 



and g 



sinh[7rTg + ] 
sinh[irT£; e ] ' 



E a - E_ = E e - E+ = El±3l. -k>0, (10) 

we can prove df/dT < and dg/dT < with the inequal- 
ity y tanh x < x tanh y when x > y > 0. Then we have 
dP/dT > O.With straightforward calculations, it is also 
proved that under the condition k ± a (nb) » c (nb), we 
have dP/dnb > [16]. Therefore, the molecular proba- 
bility P increases with T and the molecular number nb- 
Therefore, in order to obtain a high molecular probability 
P, the speed of the variation of the detuning A should 
be low enough. In Fig. 2, the behavior of P as a function 
of T is shown. It can be seen from this figure that when 
nb = and T is larger than 25\~ 1 , P is approximately 
unit and a molecule is created with unit probability from 
two atoms. However, when n b — 5, we have P w 1 when 
T is larger than 5x _1 . 

FIG. 2. The molecular probability P as a function of T in 
case the molecular number n b — 0, 2 and 5. Here, we have 
assumed k — 20\, (U x — Ub/2) — 0.5x and the unit of T is 
X -1 . It is obviously that P increases with T and n b . 



The influences of the amplitude k of the detunning 
and the atomic scattering coefficient (U x — Ub/2) on the 
molecular probability P is shown in Fig. 3. It can be seen 
that with the same operation time T and molecular num- 
ber n b , as the amplitude k increases, the molecular prob- 
ability P decreases. On the other hand, when the scat- 
tering strength is small enough so that (U x — Ub/2) j\ 
< 10, P has almost the same value with respect to dif- 
ferent values of (U x — Ub/2), i.e. the effect of atomic 
scattering can be ignored. When (U x — Ub/2) becomes 
larger, P increases with (U x — Ub/2). It is pointed out 
that in Fig. 3 we assume k ~ 10 2 and T ~ 10~ 4 so 
that e kT < 1. Therefore, the figure is drawned in the 
domain where the Landau- Zener formula is not applica- 
ble. In case of larger k or T, the same conclusion can be 
obtained with the expression of P in Eq. (9) given by 
Landau-Zener formula. 

FIG. 3. The molecular probability P as a function of k 
with different value of (U x — Ub/2). In the figure, we note 
m = (U x - U b /2) . We assume n b = 1 and T = lCT^ 1 - The 
unit of m and k are all X- The unit of P is 1CP 3 . 

Apparently, the quantum state in Eq. (7) can also be 
considered as a atom-molecule entanglement state writ- 
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ten in the Schmit decomposition form. The entanglement 
of this state can be described by the entropy of entangle- 
ment S = — J2i=i 2 \Di { n b)\ In \Di (nh)| as a function of 
the Schmit numbers \D\ (rib)\ and I-D2 (~ifc)|- The varia- 
tion of S as a function of T has been shown in Fig. 4. 
It can be seen that the entropy S has a non-zero value 
when T = 0. This is because in the limit T = 0, we have 
P = 1 — E + E_/ (E a E e ) 7^ 0. As T increases, S increases 
to its maximal value 0.49. In this case, the finial state is 
just the maximal atom-molecule entanglement sate 

\Bell) = ^=[\e,n b ) + \g,n b + l)}. (11) 

After crossing its maximal value, the entropy S slowly 
decreases to its asymptotic value in the limit T — > 00. 

FIG. 4. The entropy S of entanglement as a function 
of T with different n b . Here, we assume k = 10x and 
(U x - U b /2) = 0.5x- The unit of T is x" 1 - 

Now we are in a position to cosider the realization 
of steady coherent molecular output with the nonadia- 
batic transition method. To this end we assume a plused 
fcrmion beam with pluse lenght T and period t between 
each two pluses is used in our system. We assume in 
each inerval t the two Raman laser beams is turnned on 
and the detuning A is changed with time as in Eq. (3). 
Therefore, if the state of the total system is just \e,rib) 
before the Raman beams are turnned on, the transition 
from state \e,n>b) to the state in Eq. (7) will happen 
during the time interval r. In this case, if the loss of 
molecule is considered, it is known that the evolution of 
the molecular population of the molecular field can be 
described by the master equation [14]: 



d 

dt Pn 



1-1^2 (n)l' 



Pn 



r a \D 2 (n-l)\ 2 p n ^ (12) 



-jnp n + j{n+ l)p 



n+l- 



where p n is the propobality that there are n molecules in 
the lattice site, 7 the molecular dissipation rate and r a — 
(T + r) is the pumping rate of the fermionic atoms. 
Here, the coefficient \D 2 (n)| is defined as 

\D 2 (n)| 2 = l-|^i HI" 

where \Di (n)\ 2 is defined in Eq. (8). 
The steady state solution of Eq. (12) is 



enough (i.e. the detuning A is changed slowly enough) 
so that the state \e, rib) is adiabaticlly changed into state 
\g, rib + 1) during each time interval r. Therefore, we 
have I-D2 (n — 1)| 2 ~ 1 which do not depend on both r 
and n. In this case the molecular distribution p n in Eq. 
(13) is propotional to N" x /n\. It is apparently that p n is 
just the Possion distribution which is just the population 
distribution of a molecular coherent state. The average 
molecular number (n) is just N ex which can be controlled 
easily. Then we can obtain a coherent output of molecule 
with the nonadiabatic transition method. The coherence 
of the output molecular field can be described by the line 
width D. This line width can be calculated with the 
standard method as in the laser theory [17] and we have 
the result D ss 9-y/N ex . On the other hand, when the 
parameter T is not very large, it can be obtained with 
numerical simulation that we have (n) < N ex and the 
Mandel Q factor is larger than 1. Therefore, we can get 
superpossion distribution in this case. 

In this paper we have shown that in the photoassocia- 
tion process, we can create the atom-molecule entangle- 
ment state with steady atomic and molecular probability 
by varying the detunning of the Ramman laser beams. 
The influence of the variation parameters of the detun- 
ning on the molecular probability P have also been inves- 
tigated. As we have shown, if we want to create a high 
probability P, the operation time {AT) should be long 
enough. However, it can be seen from Fig. 4 that to ob- 
tain the maximal entanglement state in Eq. (11), the op- 
eration time needed is much shorter than the one needed 
to create a pure molecule state. We have also shown 
that steady coherent output with Possion or superpos- 
sion molecular distribution can be obtained by turning 
on and off the Raman laser beams for many times. It is 
pointed out that, since in photoassociation process the 
"Rabi frequency" x 1S much larger than the molecular 
dissipation rate 7 [14] , the condition that the time inter- 
val r < 10 2 x _1 << 7~ x can be realized. Therefore, the 
scheme in our paper can be used in the atom-molecule 
wave function engineering. 
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T^r \D* (n~l)\ 2 N ex 

Ml 1 



(13) 



1=1 



where po is the normalization constant and N ex = (r a /j). 
Since the value \D 2 (n — 1)| 2 is insensitive with the inter- 
action time t, the steady molecular distribution p n is 
also insensitive with r. A special case is that T is long 



a Electronic address: suncp@itp.ac.cn 

b Internet www site: http:// www.itp.ac.cn/~suncp 



4 



[1] u. N. Dcmkov and M. Kunike, Vestn. Leningr. Univ. Fis. 

Khim. 16, 39 (1969). 
[2] K. A. Suominen and B. M. Garraway, Phys. Rev. A 45, 

374 (1992). 

[3] P. S. Julienne et ai, Phys. Rev. A 58, R797 (1998). 
[4] D. J. Heinzen et ai, Phys. Rev. Lett. 84, 5029 (2000). 
[5] Roahn Wynar et ai, Science 287, 1016 (2000). 
[6] C. Mckenzie et ai, Phys. Rev. Lett. 88, 120403 (2002) 
[7] J. M. Gerton, et al, Nature (London) 417, 529 (2002). 
[8] C. A. Regal, et ai, Nature (London) 424, 47 (2003). 
[9] S. Durr et ai, Phys. Rev. Lett. 92, 020406 (2004). 
[10] C. P. Search and P. Meystre, cond-mat/0404006 
[11] M. O. Mewes et ai, Phys. Rev. Lett. 78 582 (1997). 
[12] L. D. Landau, Phys. Z. Sowjetunion 2, 46 (1932); C. 

Zener, Proc. R. Soc. London Ser. A 137, 696 (1932). 
[13] J. R. Rubbmark et ai, Phys. Rev. A 23, 3107 (1981). 
[14] C. P. Search et ai, Phys. Rev. Lett. 91, 190401 (2003). 
[15] P. W. Anderson, Phys. Rev. 112, 1900 (1958). 
[16] It is easy to prove that E e , E a and E± satisfy dn b /d(E e + 
E a ) > and E + dE e /dn b > E e dE + /dn b .Considering Eq. 
(10) and the inequality y tanha; < xt&nhy when x > y > 
0, it can be seen that these inequalities lead to the results 
df/dn b < and dg/dn b < 0. Then we have the result 
dP/dn b > 0. 

[17] M. O. Scully and M. S. Zubiary, Quantum Optics (Cam- 
bridge University Press, Cambridge, U.K., 1997). 



5 



